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Abstract 

The nonperturbative parton distributions, obtained from the Lorentz 
contracted wave functions, are analyzed in the formalism of many- 
particle Fock components and their properties are compared to the 
standard perturbative distributions. We show that the collinear and 
IR divergencies specific for perturbative evolution treatment are ab¬ 
sent in the nonperturbative version, however for large momentaj p? 3> 
a (string tension), the bremsstrahlung kinematics is restored. A pre¬ 
liminary discussion of possible nonperturbative effects in DIS and high 
energy scattering is given, including in particular a possible role of 
multihybrid states in creating ridge-type effects. 


1 Introduction 

The standard partonic model [D El Elia E] is a basic element of the modern 
QCD and is widely used in the treatment of high-energy processes, see mmn 
for a recent review. 

The main principle of the standard partonic model is the assumption of 
the almost free relativistic motion of quarks and gluons with only perturba¬ 
tive interaction mechanism of their evolution, given by the DGLAP equations 
PEQun]. In this process, however, as a necessary element occur IR and 
collinear singularities, which are cured by making appropriate cutoffs, imply¬ 
ing that perturbative laws are effective everywhere in the considered region. 
The neglect of nonperturbative (np) treatment is due to several reasons. First 


1 


of all, the processes under consideration are in the systems, moving with high 
velocity, and till recent time one could treat those only in the perturbative 
manner, using Feynman diagrams for free elementary constituents. Secondly, 
the study of the behavior of high-excited systems even in the c.m. system in 
the np approach is not well developed, and the standard np rules here are 
absent. 

An approach, where perturbative and np methods to the processes with 
high momentum particles were combined, namely the factorization method, 
is well developed, but needs some corrections [I21[I3]. Therefore there is a 
strong need for the development of the np side of the high energy QCD, using 
in particular methods available for relativistic np treatment. 

During last decades there had been serious efforts in formulating the path 
integral methods for QCD based on a general dynamics, including pertur¬ 
bative and np contributions [HI Ea [m [m [m [m [20]. Recently a most 
advanced version, including also external QED helds, has appeared [2T1 [22], 
and the spin-dependent dynamics was treated on the same ground [23]. 

Basing on this powerful method, the author was able recently to consider 
the QCD dynamics in the fast moving system, and has discovered, that 
exploiting the well-known Lorentz contraction effect one immediately obtains 
the parton scaled distribution from the original np wave function in the c.m. 
system [23]. In this way the knowledge of the np wave function, or more 
generally, of the total Fock column of the system in question allows to write 
down the corresponding parton distribution functions (pdf). 

Explicit examples of the unpolarized valence and sea quark and gluon pdf 
have been given in [25] in comparison with those calculated before in [6l|7l|8] 
and the similarity of the obtained pdf with the pdf from [6l[7l[8] is striking. 

Moreover it was shown in [23], that the polarized pdf of the proton can be 
obtained from the np proton wave functions calculated in the full three-quark 
Dirac bispinor formalism (the so-called factorization approach). 

The important point stressed in [23], is that pdf can be obtained from 
the Eock components both in the c.m. and in the fast moving Lorentz frame 
and the Eock composition coefficients (relative probabilities of components) 
are boost invariant. 

Moreover for the polarized pdf it was suggested in [23] to use the mul¬ 
tihybrid model for the excited DIS states, which yields ^(Es) = 0.182 for 

= 10 GeV^ which is not far from the data (respectively 0.26 for = 3 
GeV2). 

Thus it might be a possible way for the solution of the old proton spin 
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problem [221 EZl I2B]- Moreover, the DIS generated polarized pdf are argued 
to refer not to the ground state proton, but rather to a high excited baryon 
state with an indehnite spin and large gluon admixture. 

In the present paper we write down explicit connection of pdf with Fock 
components and stress the possibly important role of multihybrid states for 
pdf of high excited hadrons. We henceforth study the Hamiltonian and 
wave functions of hybrids and multihybrids and their contribution to pdf’s. 
Finally we turn to the evolution process and make the np analysis of IR 
and collinear effects, which are hnite in the np framework and divergent in 
the perturbative formalism. We demonstrate how the perturbative picture 
develops at high momenta, p ^ s/a, and how it is goes over into the np 
regime at low momenta. 

The paper is organized as follows. In section 2 the Fock component 
structure of pdf is outlined, in section 3 the hybrid and multihybrid states 
are analyzed and the gluon and quark pdf are calculated in comparison with 
known data. In section 4 the IR and collinear effects are studied both in the 
perturbative and np language. Section 5 contains summary and discussion. 

Two appendices contain details of derivation. In appendix 1 the path- 
integral calculation of the standard IR and collinear amplitudes demonstrates 
the absence of divergencies in the np dynamics. The appendix 2 contains cal¬ 
culation of the triangle and quadratic Feynman diagrams with conhnement, 
which shows how perturbative dynamics occurs for large momenta, p^ S> a. 

2 Parton distributions from Fock components 

We follow here the formalism introduced in [25l [311 [32] and write the total 
Hamiltonian and the Fock wave function as follows 
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( 1 ) 


Tw — '>Pn,{k} 


( 2 ) 


m{fc} 


Here m{k}, ^ refer to the number and types of constitients {q, q, g), modes 
of excitations, spin and momentum projections etc. 
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The orthogonality conditions look like 



/ '^N'^udV = Y = dNM- 

m{k} 

(3) 


j '^m{k}'^ri{p}dT = 5mnd{k}{p}- 

(4) 

From ([T]) 

one hnds the equation for cYk] 



^n{p}{dd^N — dd^n{p]) — Y^ ^m{kYn{p},m{k} 

m{/c} 

(5) 

where 


(6) 


The Fock component approach was exploited earlier in [29], see [30] for 
more references, in the framework of the light-cone formalism and perturba¬ 
tive dynamics for the polarized parton distributions. Below we shall handle 
the np components in the unpolarized case, and in the next section specih- 
cally for the multihybrid states. 

In what follows we shall have in mind the large Nc limit, classifying 
different contributions to namely in the leading order all states are made 
of gluons only (glueballs) or of hybrids: q, q plus any number of gluons. In 
the next order of the 1/iVc expansion the decay into two bound states due to 
the process g ^ qq is possible. 

We are using the following normalization condition as in [25] for the n 
particle bound state. 

and = Modxi, where Mq is the c.m. energy (mass) of system. Writing 
((7|) as / \(fn\‘^dTn = 1, one can write the u, d, and g distributions trough the 
Fock components (p^g, (p^qq, ^ 3 qg, as follows (4gg = uuudu) 

u{x,k\) = Un{x,k\) = \czq\^ f \<p3q\‘^dT3g Y (^i -x) 

n=3q,iqq i=l,2 

+ \CAq\^ j WAqq?dTAqqY^^‘^\^^X ~ k_L)5(Xi - x) + (8) 

i=l 
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( 9 ) 


+ \cm 9 ? / Y. - kx)5(a;i - x) + ..., 

i=l,2 

U{x,kl) = J \(fiqq\‘^dT4qqS^^\'kf’ - k±)6 {X5 - x) + ..., 

g{x, kl) == J \(p 3 qg\‘^ 6 Tsqg 6 ^^\'k^l^ - k_L)(5(a;4 - x) + .... (10) 

Assuming equal contributions of additional q and q, one can write the 
standard relations 


u'^{x,k‘]_) = u{x,k\) — u{x,k‘]_) (11) 

J u'^ix, k‘i)(fk±dx = 2 (|c 3 g|^ + |c4gg-|^ + + ...) = 2, (12) 

where we have used the total orthonormality condition ([3]), a similar relation 
for the d quark, d'"{x,k\) obtains from ([8]), ([9]), replacing the sums J 2 i=i ,2 
and by the terms with i = 3 and i = 4 respectively, resulting in the 

normalization equation 


J d^{x, k]_)d^k^dx = \C3q\^ + \CAqq\^ + \C3qg\^ + .... = 1 (13) 


From (El), (IT^ one easily obtains the Gross-Llevellin-Smith and Adler 
relations for and . 

We now turn to the momentum relations, and to this end we make evident 
the X dependence of each pdf, namely, as shown in [25], Eqs. (17), (18), the 
dependence of (p on , Xi enters in the form, where there are present the 
c.m. values of the Ath longitudinal momentum and the c.m. energy of 
the Ath quark, antiquark or gluon, namely 


„(o) , „_(o) 

x, = ^= \ , Px/r^ = Mo 


P P^/T 


or 


P\\ 


(0) 


= Mn Xi - 


M' 

H 

Mn 


, U —)■ 1 


and writing in ip only momenta of the Ath particle, one has 

- Ui)), 


(14) 

(15) 


(16) 
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where z/j = ^ is a part of the total mass contributed by the z-th c.m. energy. 

The |(pp depends on Xi in the contribution + (Mo(a:j — Vi)Y = pj- 

The parton distribution in the hadron is expressed via the np boundstate 
wave function ipisf as 


DUx, kx)= Y. n {Y kx) i (i 

X 


X Mo^|(p7vp^(5^^)(k_L - k_Lj)(5(a; - Xj) 
j 


(17) 

and it satishes the normalization conditions 



J d‘^k±dxDjj{x, k±) = 


(18) 

J d‘^k±xdxDff{x, k_L) = 1. 


(19) 


Taking into account (IT^ . (IT^ . (IT^ . one can write for the proton pdf’s 
the normalization condition 


[ dxx[u{x) + d{x) + u{x) + d{x) + g{x)] = 1. (20) 

Jo 

We now turn to the np description of the Fock components (pAr(ki^,,,k„). 


3 Multihybrid Fock components of a hadron 

It is clear, that the general equation ([5]) for the Fock components 

= E (21) 

m{k} 

contains different stages of evolution n{p}, each of these is the one or more 
bound states of q, q and g, i.e. of mesons, baryons, glueballs, and hybrids. 
In the lowest order of the expansion in l/W, initial (primordial) mesons and 
baryons are connected only to mesohybrids and baryohybrids respectively, 
containing arbitrary number of gluons. In the next order of 1/W one of the 
gluons can create the qq pair and split the total bound state into hadrons 
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Figure 1: The multihybrid meson conhguration with n gluons inside a con- 
hning string 


etc. Therefore it is of interest to study the spectrum and wave functions of 
multihybrids and we shall do it for a mesonic multihybrid, see [33l [SU [35l [36] 
for earlier discussion of one-gluon hybrids, [3zi ESI ESI snumiisiisn] for the 
theory of hybrids in the held-correlator formalism, [H] 03 SSI SZl ESI ?, 
EDESIESI for the treatment in the framework of potential and QCD sum 
rule methods, [Ml [55] for lattice data, and [56] for a review. 

We consider as in [311 ES] (Fig. 1) the system of n gluons on the string, 
connecting quark (z = 1) and antiquark {i = n + 2) and we are using as it is 
usual for the path-integral Hamiltonian method [211 ESI ES] , the einbein form, 
with effective masses Ua, oJb for qq and cu*, z = 1,... n for gluons and izfc- param¬ 
eters, replacing the linear conhnement pieces between q,g, ...,q, cr|rfc|, k = 

1,... n + 1, by the quadratic terms cr|rfc| ^ | , and Uk to be found 

from the stationary point analysis of the total energy - this procedure is 
known to yield accuracy of the order 5 10%. 


The total Hamiltonian (without gluon exchange and spin-dependent terms) 

is 


g(P)^ + 

' ' 2oJi 2 2a 

t=a,b * 


2 n +1 ^2 


P? 


-E- + E 


( 22 ) 


where a = J2i=i Ui, Vi are to be found from the stationary 

point conditions 


dH _dH dE(P) _ dE(P) 

dui dvi dui dvi 


(23) 
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It is easy to see, that using the l.h.s. of fl23p one arrives at the “relativistic 
model of hybrids” in the sense of Godfry and Isgur 0. we shall instead use 
the r.h.s. form, which will give us a simple result for -E(P) and wave functions 
(the “einbein method”). 

We start with the c.m. system, P = 0, and introduce a trial wave function 


d'o = Nexp 


n+1 


fc=i 


n+l 

]v=n 


k=l 



( 24 ) 


Leaving details of the calculation to the appendix 1, we write below the 
result of the minimization (denoted with the superscript (0)) with respect to 
{i>k,uji, CKfc} in the limit nia = mi, = 0, n ^ oo 





^ 0.62 GeV 


= 1 ^ 0 ) 



a 

275 


( 25 ) 


= 2a.„n + 0(l/n) = 1,24 n GeV. (26) 

The average intergluon distance is found to be 

= 0.68 fm. (27) 

v 2cr 

The value fl2^ is an upper limit, and the actual mass should be lower. To 
check the accuracy of our calculations we consider the case of the simplest 
hybrid state: qqg, in which case the same type of the trial wave function 
yields M^) = 2.1 GeV, while more accurate calculations [3Z] give the lowest 
mass ^ 1.54 GeV. 

Hence one can conclude, that the energy interval between the states of n 
and (n + 1) hybrids (one additional gluon) is 



AM = - M(”) ^ 0.8 = 1 GeV. (28) 

The energy intervals of the transverse vibrational modes and 

longitudinal vibrational modes have been calculated in [37] last 

reference, for the case of the gluon string with hxed ends and are equal 


A p(LVM) _ 0-53 GeV _ 0.69 GeV 

[i?(fm)]V3’ i?(fni) 


( 29 ) 


From fl2^ one can conclude, that the lowest types of excitations for long 
multigluon strings are of the TVM and LVN types, while spin-spin interaction 












interval AMgg is of the order of 0.3 GeV. It is understandable, that for large 
n there appear collective excitations of the TVM and LVM types. 

We now turn to the partonic form of our n-gluon hybrid wave function, 
from fl2^ one has in the fc-space 

^5,(k,....k„_.) = ^(-) A' = n + 2. (30) 

Here (pn is normalized as in([7]), Mq is the total energy in the c.m. system. 
From fl30l) one can deduce that the wave function for the Gth degree 
of freedom in fl30|) can be written in the form of (1161) = ^2(kl + 

and inserted in flTTD to produce the pdf of the multihybrid 
state. Integrating over (ik_i_j dxi in flTTD one hnally obtains with the notations 
= 2a, Mq —?• Mjv the form of pdf suggested and studied in |25] 

fi^\MN\x - v\) = ^J^^^exp (-^{x - Ui)A (31) 

K \ J 

where i = g,q] is dehned by normalization condition 

/ fi^\x)dx = l (32) 

Jo 

and we have introduced effective masses of quarks and gluons Mjv = 3a; + 
Nm, 3a; = Mp, so that 


I'n 


m u 

Af/v' M]^ 


(33) 


Using fl3T|) one can proceed as in [25] to calculate the pdf’s of gluons g{x) 
and valence u quarks u{x), and we shall neglect as in [25] the contribution 
of the nucleon itself, so that in multihybrids iV > 1, and concentrate on the 
region of small x, x < 0.1. One has for a sequence of N- multihybrids with 
with probabilities \Cn\‘^ 


uJx) = 2 


\T _ 1 ^ 


exp 


N=1 
oo 


g„(x) = i: |C„|yf>Jvl^^±^exp 

N=1 ^ 




4, 


- 

(34) 

- ■ 

(35) 
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In what follows we shall slightly generalize the discussion in [2S] and 
assume that \C]^\‘^ decrease as 1 < 7 < 2 , 

= \Cm? = N-^\c\\ |c|2 = C(7), (36) 

where ^( 7 ) is the Riemann dzeta-function, = (2,62)“h 

As a result one obtains for g[x) and u{x) 

g{x) = \c\^{aix^~‘^ + a 2 X^~‘^), (37) 

Un{x) = \c\^{bix^~^ + b2x'^~'^ + 60 ), (38) 

where ai,bi depend on parameters K,m{uj = ^ = 0.31(GeV). We can hx 
these parameters, as in [25] , using the explicit form of the lowest hybrid with 
= 1 from [321 EZ], K, = Y ~ 0.313 GeV. 

In this way we have the only htting parameter 7 , which we can choose 
to correspond roughly to the gluon pdf at = (3 10) GeV^. In this 

way we can consider g{x) and u{x) as starting pdf in the DGLAP evolution 
equations instead of standard initial state functions with roughly 15-20 free 
parameters. Following [25] we take 7 = 3/2 and obtain 

ai = 1.5, 02 = 1, bi= 2.82, &2 = 2.70. (39) 

As a result one obtains xg{x) at the points x = 10 “^; 10 “^ to be equal to 
3.87; 1.38, which is close to the PDG data at = 10 GeV^, respectively 5; 
1.5 (cf. Fig 19.4 in the last ref. of [6l[7l|8]). 

The same type of similarity between fl38|) and and PDG data exists for 
xUy{x), which implies, that the sequence of multihybrid states can be used 
as a part of evolution process. It is also interesting that g{x) contains an 
additional power of x~^ as compared to Uy{x), which also nicely corresponds 
to the data. One can now check whether multihybrids can produce reasonable 
amount of sea quark pdf. This check was done in [25], where it was shown 
that the qq DGLAP evolution of fl37|) in the interval (1-10) GeV^ produces 
u{x, 10 GeV^) 0.04g{x, 10 GeV^), in good agreement with PDG data [ 6 ]. 

As a hnal check we compute g{x,Q'^ = 10"^ GeV^), as a result of the 
DGLAP evolution of our form fl37|) . fl59D for = 10 GeV^ and obtain 

= ^^(p(a:), (p(x) ~ |cp(0.917x“^/^ -|- x~^ — 3Ax~^^‘^) (40) 

alnq)"’ tt 


10 




which yields roughly xg{x = 0.01,10"^ GeV^) ~ 3.5, whereas the PDG data 
gives a 2.5 larger value. This might point out to an additional mechanism 
(e.g. of the BFKL type) at very large and small x. 

We have considered above only the region of small x, x < 0.1. Larger x 
need some modifications. Namely, introducing the quark-counting rule fac¬ 
tor (1 — xY^^'^ into (Ell), where p{N) is growing with N, one obtains 
modified normalization factors and an effective upper limit in the sum¬ 
mation over N, /Gnax - X, since M%{eff) s. These replacements 

do not essentially change our results for iVmax ^ 1 , but are essential for very 
small X, where the contribution of the multihybrid contribution is strongly 
damped by p{N) as compared to the nongluonic excited nucleon states. 

Goncluding this section, we have shown that the multihybrid mechanism 
at the initial or intermediate stage of parton evolution can be considered on 
the same ground as other realistic contributions and may have its own part 
in the final reaction products. 

Note, that our small x asymptotics, Ug{x) ~ x~'^, is different from the 
reggeon-type asymptotics, and is a result of the multihybrid state formation 
in the general multigluon ladder-type diagrams. It is possible, that the latter 
provide both types of asymptotics, the more common “hard” momentum 
distribution when each loop in the ladder with a large rapidity ratio yields 
(^aln^^, and a soft (or coherent) momentum distribution, as in fl30l) . which 
possibly occurs for very large s. We study below in the next sections and 
appendix these two regimes in the example of the square box diagram with 
confinement, and discover two possible regimes, depending on whether the 
external momenta are much larger, or comparable with the string tension. 

Significantly as one of possible signals one can consider the ridge-type 
configurations [721 [731 [71], where experimental signals are along a straight 
line, which mimicks the average form of a large multihybrid with iV S> 1. 


4 Nonperturbative vs perturbative regimes 

The purpose of the present section (and the Appendix 2 ) is to demonstrate, 
that IR and collinear singularities are absent in QGD (in contrast to QED), 
which however does not imply the absence of the logarithmic and double 
logarithmic terms in the perturbative series. 

As a standard approach to the evolution and cross sections in high-energy 
QGD one is using the QED processes of e ^ e'y, 7 e~^e~ with the IR and 
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collinear singularities and perform the analogous perturbative procedure in 
QCD, leading to the singularities e.g. in fragmentation functions and average 
number of partons etc. However, in QCD both IR and collinear singularities 
are absent, as we demonstrate below, and the multiplicative gluon emission 
leads to the formation of heavy multihybrid states, in contrast to the QED, 
where the multiphoton emission is a real background process, which is treated 
introducing an experimental bound - the minimal energy Ei of the registered 
photons [58l Ell [60] . 

As a consequence in QED the observed cross section of the high process 


contains Sudakov double-logarithmic corrections of the type ^ In 



In QCD the minimal energy Ei is absent, since gluons are never emitted 
as free particles and its emission inside hybrid states costs an increase in the 
total mass of ground states AE 0(1 GeV). Moreover, the quark mass niq 
cannot enter in the hnal equations due to conhnement for niq < \fa = 0.42 
GeV, hence the multigluon emission might have another structure. In the 
DGLAP evolution equations this difficulty is avoided considering evolution 
of quark, antiquark and gluon pdf’s ff{x, Q) as functions of momentum scale 
Q, with the kernels ensuring logarithmic dependence on (a weak violation 
of the Bjorken scaling), and on x at small x. 

As we argue below in this section, and in the Appendix 2, in QGD the 
IR and collinear singularities are absent, but the perturbative series terms 
can indeed have the known logarithmic form, and the role of the cut-off 
parameter is played by the string tension, so that at high momenta, S> a, 
the perturbative kinematics prevails. 

We start with the standard perturbative picture. To put the problem 
in the most simple form, consider the process e+e“ —)■ qqg, with the cross 
section for free q, q,g equal to (see e.g. |5], chapter 17) 


d(T , 
dxidx2 


mg) = 



2 as x\ + x\ 

Svr (1 - a;i)(l - X 2 ) 


(41) 


where xi, 0 : 2 , 0:3 are ratios of q^q and g energies to the e'''e energy ^/s and 
one can write 


E E 

1 - xi = a; 2 ^(l - cos 6 ' 2 g), 1 - 0:2 = a;i^(l - cos 6 'ig), (42) 
v s v ^ 


and these denominators arising from the virtual propagators of quarks: 

(p'l)^ — ml = (pi +Pg)‘^ = m‘l + s/sEgXi^l — cosi^). One can see the collinear 
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Figure 2: The Feynman diagrams for the qq Green’s function, which con¬ 
tribute to the qqg cross production cross section 

and IR singularities of the unconhned QCD, which are essential element in 
the balance of the integrals dehning the terms of the perturbative series. 


In general, the bremsstrahlung probability can be written as 


dPf 


brems — 


aC^ cPk± dx 

TT^ k\ X ^ 


(43) 


where x ^ —, or else x ^ ^ for the deep inelastic evolution. In the series of 
the gluon cascade amplitude (the Regge-Gribov amplitude) one has for the 
term with n internal lines 


dxn P dxn-i P dxi 1 1 

- -(ca^ln-) 


/ \n f dXn f dXn-l f 

/ — / TT—- / 


1X2 Xl 


ni 


X 


which yields for the total sum the standard answer 


^ 4 ~ exp (ca^ In -) ~ 

n ^ 



( 44 ) 


( 45 ) 


One can see, that the crucial property of the Regge-type behavior (I45|) 
is the bremsstrahlung-type energy distribution at each step of the gluon 
cascade, as it is also in fl4ip . Therefore it is interesting to calculate the same 
cross section fITT]) with conhnement taken into account. 

To make a comparison with the np (conhned) QGD, one must calculate 
the diagrams, shown in Fig.2, and one has in mind, that all the area inside 
the outer contour of the diagrams, is covered with the conhning him (shown 
by the vertical lines in Fig. 3). 
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Figure 3: The path integral representation of the qq Green’s function. The 
regions denoted by 1,2,3 refer to qq,qqg and qq intermediate states respec¬ 
tively 


We calculate the qq Green’s function in the path integral formalism in¬ 
corporating conhnement and gluon exchange in appendix 1 and write down 
the result in the form, appropriate for the Fock component language, namely 




{y) 

{2uy 


(p(0) 


F, 


12 


E — El ^‘ZoOg E 


1 F23 1 

— E2 y/Zujg E — E3 


fp{0) 


(46) 


where the gluon creation matrix elements Ffc are expressed via overlap matrix 
elements of meson and hybrid wave functions [HU |32l EU [62], and u, Ug are 
average energies (from the Hamiltonian minimization) of quarks and gluon. 
From [32] one has 


Vu 

y/IEEg 




23 




= 9 


O.OSGeV 


(47) 


To obtain the cross section e'''e —?• hybrid{qqg), one takes the disconti¬ 
nuity of fi46|) in energy due to the factor and obtains 


da{e^e —)■ hybrid) ~ 


(p(0) 


F 1 1^ 

^yZujg E2 - Ei_ 


6{E - E2) 


(48) 


where F ^2 — ~ 0.8 1 GeV for Ei ^ rup, E 2 = Mhybrid and one does not 

have any IR or collinear singularities in agreement with our discussion and 
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the beginning of this section. It is also essential, that the strnctnre of the 
hnal states in 0481) in the np regime is qnite different from the pertnrbative 
QCD, namely for low energy resolution one should sum up over all collective 
states of the hybrid, and in the case of high energy and multihybrid state 
the overall entropy of states can be very large. 

To resolve the apparent contradiction one can consider the Feynman di¬ 
agrams, shown in Figs. 4, 5, where we take into account conhnement as the 
area law of the Wilson loop, e.g. formed by the 4 sides of the rectangular in 
Fig. 5. 


W{C) = exp(-crS'™,„(l, 2 , 3,4)) (49) 

As shown in the appendix 2 , the approximate answer for this diagram is 

GM = n / dUi r 6«) (50) 

*=i -^0 

and 

=Pi-qi + qk (51) 

where qi, q^ are momenta on the adjacent lines in the vertex i. The property 
of Jo- is such, that at large Pi,qk, (5*-*^)^ 3> a, it becomes a product of the 
vertex 6 - functions 


Jo((6(*))2 >a) ~ (52) 

i=l 

As a result G^i^pi) coincides in this limit with the standard perturbative 
answer, where the bremsstrahlung property flT3l) appears, but now for large 
(as compared with y/a) momenta. 

At the same time for low momenta, ( 6 ^*^)^ < a, one has four d^q integra¬ 
tions, which exclude any possibility of IR or collinear singularities. 

Thus one may imagine for a sequence of ladder-type diagrams two types 
of asymptotics, depending on the internal momenta distributions: a) (IT7|) . 
(I48|) for perturbatively ordered momenta; b)multihybrid-type asymptotics 
considered in the previous section for np ordered momenta, |APj| ~ y/a. 

We are now coming to a possible general picture of np evolution in DIS 
or high energy hadron reactions. We assume, that in the intermediate or 
initial state there appears a multihybrid state (or a family of multihybrid 
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states with close-by masses). This multihybrid state originally (by evolu¬ 
tion) is perturbatively (bremsstrahlung-type) ordered, which corresponds to 
the reggeon-like asymptotics and the BFKL-type desctiption, but has a pos¬ 
sibility to develop into a quark stable multihybrid state, described in the 
previous section. It subsequently decays into more multihybrids and hnally 
into qq and 3g hadrons, due to perturbative and np string breaking. In this 
way the gluons contribute the most amount of the collision energy at this 
primordial stage, which in the evolution process is transferred to the hnal 
hadrons. 


5 Summary and discussion 

We have tried above to formulate main features of the np approach to the high 
energy production of hadrons in DIS or hadron-hadron collisions, which might 
be an alternative or a complement to the standard QCD picture 
The basis of our approach is the fact, that the Lorentz contraction rule of 
wave functions automatically leads to the parton-like form of its momentum 
depenence [2lj, and one the c.m. wave function to rewrite it in terms of 
partonic variables, p_|_ and x, namely = pi + Mq[x — z/)^, where Mq is 
the c.m. energy of the object. 

This allows to write down the pdf’s of the object, e.g. of the fast proton in 
terms of the c.m. wave function, as it was done in |2S] both in the polarized 
and unpolarized cases. In this way one obtains the np components of the 
partonic set of pdf’s, and the point is how to represent wave function of the 
arbitrary moving complex object and subsequently how to extract from it 
the total set of pdf’s: valence and sea for any quark flavor. 

This is done at the beginning of the present paper with the help of Fock 
components of the Fock column wave function, where each line corresponds 
to a dehnite ensemble of g, g, g with the same total quantum numbers. We 
have expressed pdf’s as a sum over Fock components, satisfying usual normal¬ 
ization conditions. We have also shown before in [25], that the normalization 
of the Fock components is boost invariant, and hence it is possible and often 
convenient to calculate pdf’s from the Fock column in the rest frame. 

Using that, we have considered in section 3 the basic building block of 
the high excited Fock column - the multihybrid state, which may be one of 
the basic primordial state of the high energy evolution of the Fock column 
in the 1/W expansion, where the g ^ qq decays are 0{1/Nc), as compared 


16 


to gluon creation and absorption. 

We have found the wave function of the multihybrid and discovered that 
it is highly correlated with an average energy per quark and gluon around 
0.6 GeV and average intergluon momentum |Apj^| around 0.36 GeV, and 
\Axi\ ^ 5jV ^ where N is the total number of particles. 

Such a complicated object for large N has a very large enthropy and can 
be stable enough with respect to emission and reabsorption of quarks and 
gluons, and hence can be a prototype of the structures of the ridge type, 
observed in pp and AA collisions [651166], [67]. 

We have also shown, that the appearance of a multihybrid at some stage 
of the evolution does not violate the observed pdf’s of gluons, valence quarks 
and sea quarks, however at very large ~ Q‘^/x one possibly needs com¬ 
bined effects of both reggeons and multihybrids. 

Moreover, an ensemble of multihybrids yielding a balance of perturbative 
and np features different from the standard approach [681 EH], can be a rea¬ 
sonable alternative saturation state of the system evolution ( see [70] for a 
discussion), since for a given excitation energy E of the system the set of 
equations for Fock components, leads to the effective cut-off of the number 
N of constituents, N ^ j- ^ os^GeV subsequently hadronized. 

Finally, one should stress some similarity between the multihybrid mecha¬ 
nism and the Lund string model EJ. We note also, that the decay of a 
multihybrid can proceed or via gluon decay g qq, or else via np decay 
of the string piece into qq, which initiate different genealogical chains. Note 
also, that the strong np decay into qq pair can be accompanied by pion emis¬ 
sion, - the phenomenon observed e.g. in the heavy quarkonia [7211731171], 
which can produce pions at the earliest stages of evolution. 

The author is grateful for useful discussions to K.G.Boreskov, B.L.Ioffe, 
O.V.Kancheli and members of the ITEP theory seminar and to I.M.Dremin 
and members of the FIAN theory seminar. The author gratefully acknowl¬ 
edges the help of I.V.Musatov in working out the material of appendix 2. The 
hnancial support of the RFBR grant 1402-00395 is gratefully acknowledged. 


Appendix 1 

Calculation of the qqg Green’s function in the path integral formalism 
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We start with the path integral form of the free Green’s function 


9{.x,y) 




°° du 


ui 


3/2 


(D^z) 


xy 


-K(co) 


where T = X 4 — and 



(Al.l) 





+ m? 

2bj 



(A1.2) 


In a similar way for the product of two spinor g, q Green’s function one 
has 


\{ml-Dl){ml-Dl))8n Jo Jo ujT 
X (Al.3) 

where 

4y = tr r(mi — Di)r{m2 — D2)', W = trexpig j A^dz^ (Al.4) 

and we have omitted for simplicity the spin-dependent terms in W. To 
implement the gluon lines in the total amplitude, we can write the path- 
integral form 


{A^{u)A^{v)) = 5^y j dse ^{D'^z)uv^{u,v) (Al.5) 

where ^{u,v) = Pexp{ig Axdz\). 

We also take into account the identities [HI [HI HSl CZl ESI ESI ED], last 
reference, 

(D^z),y = iD^zUd%iD^z)^y (A1.6) 


POO PS PT\ pOO pOO pOO 

/ ds dri dT 2 f{s,Ti,T 2 )= ds dn dT 2 f{s+Ti+T 2 ,Ti+T 2 ,T 2 ). 
Jo Jo Jo Jo Jo Jo 

(A1.7) 
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Our purpose is to write the path integral representation of the diagram 
in Fig.l. Using flAl.Sp . flAl.hp and flA1.7p one can write 

Ggq{x,y) = / dsidTidaids2dT2da2{D'^z)xud^u2ui2iD'^z)uvd^v2u23X 


x(D‘z)^,fW„e-’'^'-'‘'{D‘z')r^.cPu' 2 w[ 2 {D‘z')„,,,ifv' 2 uJ 23 {D‘z),,, Jds' e-‘^HD*z)^ 3 WA. 

(A1.8) 

One can now use the correspondence in flAl.ip to introduce the Hamilto¬ 
nians Hi, H 2 , H 3 for three sectors 1,2,3 shown in Fig.l, namely 


/ 


dsi{D'^z)xue 




e-Hdt 


u) 


(A1.9) 


and similarly for H 2 , H^. As a result the total qq Green’s function after 
integrating out the c.m. coordinates is 


Gp(T) = /G„-(x,i,)e‘P(^~y>d^(x -y) = 


_{Y) f t dUidU 2 \^3 


4: J 271 {uiU2y^^ 


e ^ P«2cni22cn 


, t' — t du[d(jj 2 


12 


Svr 


X 




If — t du'n 




Svr 


X 


T-f f du'Idu'' 


{Pv\e 


-HiiT-t') 


\ 0 )dtdt' 


(Al.lO) 


Stt J 

and hnally, integrating over Euclidean time intervals, one obtains 

Fl 2 1 V 23 1 


Ggg{E) = J Gp{T)E^^dT = 




E — El \/2ujg E — E 2 \/2ujg E — E^ 

(Al.ll) 


972 ( 0 ) 


where Vik are dehned as (see also 


Vik = = g J (pM{r) (Al.12) 


where and ^' 0 n(i'i 2 , r 23 ) are wave functions of the meson in the sector 

1 or 3, and respectively of the hybrid with /i - the gluon spin orientation. 

The form flAl.llD is exactly what is expected in the Fock component 
formalism of section 2. 
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Appendix 2 

Integral representation for the 3 and f point functions 


In this Appendix we shall consider the path integral technic [211 [22] for 
3 and 4-point amplitudes with one closed quark contour and any number of 
perturbative and nonperturbative gluon interactions. The basic technic was 
described in pil \T2\ . and we start with the Euclidean space-time. Then the 
n-point amplitude defined as 

G{Pi^\ --Pn^) = {JliPl)-JniPn)), Ji{x) = 

can be written using [21] as 


(A2.1) 




i=l 


(A2.2) 


Here we have denoted 


^'^) = PAPpewiig [ Af,dZf,)exp{g [ dria^^F^^), (A2.3) 

Jx(i-l) JO 


afii/ — 


and 


n ‘‘fff = (A2,4) 

Eq. flA2.2l) is the exact expression in QCD, when no internal quark loops are 
present, so it should be exact in the large Nc limit. 

All phase factors flA2.3ll combine in flA2.2l) into a Wilson loop factor with 
insertion of operators aP, which we denote as (hFo-), 


{W„) = (11 'I>?>.4. 


(A2.5) 


2=1 


In what follows we disregard the spin factors aF, in the first approximation 
since as it was shown in [H nsi nsi [m nHiEni ED], they give nondominant 
contribution to the amplitude at large momenta. Moreover the [mi — Di) 
factors when acting on path integrals in flA2.2j) are shown to be written as 
[mi — iqi), where g* are Minkowskian momenta on the given line (see last 
reference in [HI [TDl ESI EH ESI ESI ED] for a derivation). 
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We start here to calculate (Wj) for the triangle diagram of Fig.5, making 
the approximation of large area, Smin ^ so that 

{{W^)a ~ eyi]i{-aS^in) (A2.6) 

Here Smin = Smin{x^^\ and we assume for simplicity the straight- 

line (eikonal) geometry for the quark lines, so that the Nambu-Goto expres¬ 
sion for Smin can written as 

Smin = J^ dt dp^w^{w'y - {ww'y, lii = ^,w' = ^ (A2.7) 

and quark trajectories (from a;F) to and (from x‘d'> to are 

given as 

+ (a^® - (A2.8) 

with 

w^{l3,t) = Zf,(t)/3 + z^{t){l - 13). (A2.9) 
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Since the Nambu-Goto form flA2.7p is difficult to handle, one introduces 
as usual einbein variables and 'r]{l3,t) to write 

^-aSmin — J I ’^[w^+{^)'^w''^-2ri{ww')+ri'^(w')^]dpdt^ (A2.10) 


The form in the exponent in flA2.10p is quadratic in coordinates and 
can be easily computed to be 


— (rSrr 


V 


= / T)z/T)pexp{-a-/}, 


where 


(A2.11) 

(A2.12) 


/ = 2/i + — 2/2 - + h 2/2> 

and we have used notation v = vja and 

yx = Xyi = 1/2 = = x*-^^ — (A2.13) 

The full integral over x^*^ which is contained in flA2.2p can be written as 

/3= / 


X g*gij/i+* 92(^2 yi) *932/2 


(A2.14) 


The coordinate part of the integration in flA2.14p can be written as 

J3(a,b) = [ (A2.15) 

J Vdeta/ 

where notations are used 

an = ^; bi = -(pd) + p(3)) - + g2 

022 = |(^ + b2 = - q2 + qs 

0-12 = 021 = (A2.16) 
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Pa 


Inserting flA2.15l) into flA2.14l) one has 

h = (27r)^5(^p«) J Du J Dr^ 


d^) 


(A2.17) 


Since z/, rj are nondynamical variables, the integrals over them can be taken 
using the stationary point analysis, which amounts to solving equations 


(Jd's I _ ^ ^^3 I 

177=770 f- \d=VQ U, 


5r] 


with 


5u 

1 


4/3 = —2 In det a + —biUj^^bk- 
4(T 


(A2.18) 


(A2.19) 


Inserting u = z^o(p*'*\ <ii)^ V = Vo{p^^\ Qi) back into flA2.17p one has g*) 

and the 3-point amplitude acquires the form 

^ /* /*oo 

C!(p(i)^p(2)^p(3) =trY[ / d^qiTi(mi - iqi) / g*), 

i=i J D 

(A2.20) 

iP = (27r)'‘5(^pW)(167r)2exp(--Y/6f6f^^^M^). (A2.21) 
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The exponential factor in displays the role of conhnement on the 3- 
point Green’s function. Note hrst of all, that 61,62 dehned in flA2.16p and 
63 = + Qi — Q3 satisfy the relation 61 + 62 + 63 = 0 and therefore one can 

check that the factor in the exponent in flA2.20p is symmetric with respect 
to the replacements 61, 62 —>■ 61, 63 —>■ 62, 63; 


^(61,62) = sjbiq - (6162)2 = ^(61,63) = 5(62,63). (A 2 . 22 ) 

Secondly, the magnitude of 6 * is the measure of ’’the off-shellness” in the 
corresponding vertex. 

In the absence of conhnement 73 ( 0 , 6 ), Eq. flA2.15p goes over into the 
momentum conservation at each vertex, 

Ma,b)l^o = { 2 nfS^^\b,)S^^\b 2 ). (A 2 . 23 ) 


For nonzero a the role of a smeared 6 -functions is played by the expo¬ 
nential factor in flA2.2ip . It is clear that for large 3> ex, this factor 

strongly cuts off all conhgurations of momenta unless 6 ^ < a, z = 1,2. This 
means that in the presence of conhnement the momentum conservation at 
each vertex is fulhlled only with the absolute accuracy of the order of a, i.e. 
Qi - Qj = + 0 (\/^). 

As a result the Feynman rules for the triangle graph with the area law 
change in such a way that instead of one momentum integration inside the 
loop one has to do all three integrations over virtual momenta qi,i = 1,2,3 
with the cut-oh factor flA2.2ip . This fact strongly changes the IR properties, 
since IR divergencies are now impossible by a simple power counting. 

On the other hand, for large momenta, S> cr, qf S> ex the exponen¬ 

tial factor in flA2.2ip strongly cuts oh the diagram magnitude unless only 
one momentum direction is available, while the perpendicular directions are 
controlled by a. Consider now the 4-point function. According to the Fig.5, 
we assign the quark and antiquark trajectories as follows: 


Z^it) = X2^ + ( 0:3 - X2)^t 

Z^,{t) = Xif, + ( 0:4 - Xi)f,t; 0 <t<l (A2.24) 

and for the standard Nambu-Goto expression flA2.7ll we dehne the straight- 
line metric paths w^: 


w^{/3,t) = z^{t)/3 + - (3) (A 2 . 25 ) 
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we also denote three independent distances: 


yi = Xi2 = Xi - X2] y2 = X23 = X2 - Xs] ys = X3 - x^, y4 = x^i = X4 - Xi] 

(A2.26) 

with relation for yp 


J2yi = ^- 


(A2.27) 


i=l 


Using flA2.10p we dehne the same integral as in flA2.14D and the integral 
equivalent to flA2.15p appears to be: 

J4(a, b) = J d‘^yid‘^y2dS4exp{-[aikyiyk + hyi]}. (A2.28) 

It is convenient to extract string tension from aik, and dehne: 

u = av, ttik = crttik (A2.29) 


where 

1,1 2n- X 1,1 yx 

«11 = o “ + ^ ^ ; «22 = ^ + X - + h ^ 

2 V b b V 4 

nl.l 2 \ yx _ 1,1 2 , 

= « + «(- + ^ «i2 = tI- + h ^ - yx)] 

b b V 4 A V 

1/1 2 \- X 1,1 2\ 

ai 4 = -[- + y V- T]v); 024 = - + -(- + h X); 

A V 12 b V 

dik = dki] bi = -qi - q3-pi - Pa] 

b 2 = Q 2 — Q3 + P 3 ] b^ = q 3 — q4 + Pa- (A2.30) 

The integral J 4 in flA2.28p is Gaussian and the result is 

2 

JA{a, b) = ^^2.31) 

(det ay 


In this case as also for the 3-point function, both rj and u are to be found 
from equations equivalent to flA2.18l) . flA2.19D . 

One can simplify the foregoing expressions using the approximation, when 
the area of the (generally speaking) nonplanar quadrangle is replaced by the 
sum of areas of two triangles formed by 4 sides and one diagonal inside the 
quadrangle. This approximation is exact for the planar quadrangle and gives 
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a larger area for a nonplanar case. It can be used to obtain a qualitative 
estimate and the upper bound for the area law correction to the 4-point 
amplitude. 

With the notations from flA2.26j) one can thus write for the area of the 
quadrangle, 

sil(l, 2,3,4) « sil(l, 2,3) + sil(l, 3,4) = 


- {yw2? + \^Jvhi - [ywA?- 


The equivalent of Eq. flA2.14l) reads 

h = (27r)"5(Sp«) n 

i=l 

with 


3(4) 


'(2vr) 


(A2.32) 


(A2.33) 


bi = qi-P2-P3 + V, b2 = q2-P3 + 'P, 

^ 3=^3 + P, &4 = ^4 + P 4 + P- 

Using flA2.32j) one obtains 


(A2.34) 


h = 


(27r)^5(Sp«) e"- V^i'>l-(fe3b4)2 ^ (27r)^5(Sp®)J®- 


(A2.35) 

The hnal expression for the 4-point function G 4 (pi,P 2 ,P 3 ,P 4 ) (where the 
factor {27i)^6{J2Pi) is separated out) 


Gi{pi) = trJl [ d qiT i{mi - iQi) [ e exp{-Q). 

^=l J -^0 

(A2.36) 

The factor exp(—Q) in flA2.36l) is equal 


exp(—Q) = exp 


g^C2 

Stt^ 


'y ( ^ji qi) qj) 


i<j 


(A2.37) 


The vertex parts in flA2.37p may contain double logarithmic parts under 
the same conditions as for the 3-point function. One can see, that also in the 
4-point case one arrives at the asymptotics at large 6^ 3> a 


4=1,2,4 


(A2.38) 
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In this case the four d'^Qi integrations in (IA2.36D reduce to a single one, 
as it should be for the standard Feynman diagram without conhnement. 
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